Possible Vortex Fluid to Supersolid Transition in Solid ^He below ^75 mK 
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A detailed torsional oscillator(TO) study on a stable solid ^He sample at 49 bar with To^^O.SK, 
is reported to T below the dissipation peak at Tp. We find both the shift of period and dissipation 
hysteretic behavior below Tc ^^75 mK, with changes of AC excitation amplitude Vac- The derived 
difference of non-Unear rotational susceptibility ANLRS{T)hys across the hysteresis loop under 
systematic conditions is analyzed as a function of Vac and T. We propose that ANLRS{T)hys is 
the non-classical rotational inertia fraction, NCRIF itself, and it is actually the supersolid density 
pss of the 3D supersolid state below Tc. pss changes linearly with T down to~60 mK and then 
increases much more steeply, approaching a finite value towards T=0. We find a characteristic AC 
velocity '^40/im/s beyond which the hysteresis starts at T < Tc and a "critical AC velocity" 14 '^10 
mm/s, above which pss is completely destroyed. We obtain and Vc—h./ {rriA'io .7r)=~6-12 mm/s. 

PACS numbers: 67.80.bd, 67.25.dk, 67.25.dt, 67.85.De. 



The supersolid state, which can be characterized as a 
solid with a lattice structure yet simultaneously having 
superfluid properties, has been one of the most inter- 
esting topics in condensed matter physics. Discussions 
started in the 1960's (see recent reviews (mil) and the 
experimental search for such a state in quantum solids 
followed for 40 years till now. Unusual properties of solid 
^He have been reported for some time. The first convinc- 
ing claim was made by Kim and Chan|3|], in a report of 
non-classical rotational inertia {NCRI)^ a phenomenon 
predicted long ago by Leggett[4]. This experimental ob- 
servation with a torsional oscillator(TO) has been con- 
firmed by other groups including the present authors' 
group [HI, 0, 05 B @]- There were, however, a few funda- 
mental problems in the identification of NCRI in the re- 
ported experimental results by that time. The problems 
were discovered following newer observations, namely, 
the reported onset temperature To of 0.1- 0.5 K [9] is too 
high for the known number density of imperfections in 
solid ^He samples [10] to expect Bose-Einstein Conden- 
sation (BEC) on one hand, while the reported amount 
of the NCRI fraction goes up to ~ 20 % of the total 
^He mass [7] in conflict with Leggett's original ideas about 
the origin of supersolidity as BEC in the imperfections. 
Furthermore, the failure to observe superfiowflll, [l2[ is 
suggestive of some other phenomenon, which is actually 
responsible. Anderson proposed a picture of a vortex 
fiuid(VF) state above a real Tc [13, J.4] for the reported 
observations for the solid ^He. He took the absence of 
superflow into account and the sensitivity to the amount 
of disorder [7] and proposed the reported TO response is 
the "nonlinear rotational susceptibility", NLRS of the 
VF state, instead of NCRI for the real supersolid state 
and the latter should exist at some lower temperature be- 
low Tc. Actually another independent picture had been 
proposed by Shevchenko long ago[l5|, [3, considering a 
ID dislocation core as an origin of superfluidity. It de- 
scribes a real superfluid state in solid He below some Tc 



and still they expect something like superfluid response 
above Tc because of the dynamic properties. We do not 
consider this picture now because we do not find a direct 
connection with the experimentally found phase transi- 
tion. 

It is well-known that the VF is a state without super- 
flow, or without 3D macroscopic coherence in the case 
of underdoped(UD) cuprates[17]. It is also argued to be 
characterized by almost constant amplitude of the macro- 
scopic wave function through the real 3D macroscopic 
super-conducting(-fluid) transition at Tc, but phase fluc- 
tuations break down the macroscopic coherence [3, [3 
at Tc, which can be at much lower T than the onset tem- 
perature. To, where quantized vortices (probably in a low 
D subsystem) start to appear thermally. As a result of 
increased low D coherence length ^jj (which increases to- 
wards T= OK), a real 3D supersolid transition should oc- 
cur at some low T. A sharp peak in specific heat, which 
would indicate a 3D real phase transition, has been re- 
ported around 75 mK[19]. 

On the other hand, there also have been some attempts 
to explain features of observations in terms of classical 
dislocation motion trapped by ^He in connection to the 
response of solid ^He to shear motion apart from super- 
solid properties [20]. Actually, a similar shear modulus 
increase below about 200 mK has been observed not only 
in hep solid ^He, but also in hep solid ^He. But a TO 
response anomaly, NLRS or NCRI has been observed 
only in hep ^He[21]. 

In a recent publication we have shown [9] experimen- 
tal evidence which supports the VF picture in solid ^He 
samples at 32 bar as well as 49 bar below a common 
To ~ 0.5 K. To was determined for the first time for 
solid ^He from a detailed study of the AC excitation Vac 
dependence change at this temperature We argued 
for Vac and T dependent responses of the pre-existing 
thermal fluctuations of the phase of a mesoscopic scale 
wave function in the VF state. Stronger AC excita- 
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FIG. 1: NLRS{T) at Vac ^ 0, is displayed as a func- 
tion of The solid line through the data points is 
the Langevin function f(x)=a[{exp(bx)+exp(-bx)}/{exp(bx)- 
exp(-bx)} - l/(bx)] with a = 0.0878 ± 0.0011, and b = 0.0148 
db 0.0004. Inset shows the Vac dependence for data at each 
T^300 mK and we can safely extrapolate to Vac 0. 
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FIG. 2: The logVac dependence of NLRS of sohd ^i7e sample 
at 49 bar at constant T' s as given in the figure, obtained from 
the measurement of period change of TO. 



tions, which would cause formation of straight vortex 
lines, could suppress these thermal fluctuations of the VF 
state. Such suppression appears not only in NLRS, but 
also in the energy dissipation^. The unique log{Vac) 
linearly dependent suppression of the NLRS and its 
unique temperature dependence suggest that what we 
observed was not NCRI, expected for the 3D super- 
solid, but what is being proposed by Anderson [22| as 
NLRS of the VF state. All the observations, especially 
above the energy dissipation peak, can be well described 
by the properties of the VF stated, HI, Q. The VF 
state may have some features in common with super- 
fluid turbulence, where the "polarization" of the vortex 
tangle under rotation may be described as an ensemble 
of vortex loops and its " polarizability" under rotation 
could show similar behavior to that represented in the 
Langevin function [23^. Similar suppression of fluctua- 
tions by an external magnetic field has been reported for 
the VF state in layered superconductors [24], and cuprate 
high Tc superconductors [17]. 3D superfluidity and 3D 
vortex lines are realized in a series of 3D connected He 
monolayer systems, where 3D connectivity of the super- 
fluid is provided by the 3D connected surface of a porous 
substrate [ill. The critical velocity of this system for de- 
struction of the 3D superfluidity seems to be character- 
ized by Vc^ h/maTT, above which 3D superfluidity is de- 
stroyed and 2D superfluid features appear [i^, where a is 
the 3D vortex core size, or the minimum size of the 3D 
superfluid. 

In the present paper we describe further a detailed TO 
study on the same stable solid ^He sample at 49 bar [9], 
but extended to T lower than the dissipation peak at Tp. 
The first observation at the lower T is that the logVac 
linear slope's T dependence deviates from the reported 
l/T^ dependence at higher T for the VF state To 



make the physical meaning clearer we consider the value 
of NLRS extrapolated to Vac=0. It would represent the 
full amount of NLRS of the VF state. Fig. [1] indicates 
NLRS{T) for Vac extrapolated to as a function of 1/T^. 
It also deviates from the 1/T^ linear dependence at lower 
T and was found to follow the Langevin function quite 
well with X = 1/T^ and it approaches a finite value to- 
wards T=0, 0.088% of the total solid ^He mass. All 
the data in Fig. 1 inset were measured under a certain 
procedure, where the Vac is set to a maximum value at 
T >~0.5K and cooled down to desired T and Vac was 
swept downwards stepwise after equilibrium at each Vac- 
After completion of the measurement at a T another set is 
repeated for another T. We call this procedure "measure- 
ment under equilibrium condition". The T dependence 
may indicate that the VF state is best described by be- 
havior of an ensemble of vortex loops as was discussed 
for superfluid turbulence in [23]. Loops can be polarized 
in a manner similar to dipole systems, where instead of 
l/T^, 1/T appears. 

In this lower T range where deviation from the 1/T^ 
linear dependence occurs, we find a new feature, namely 
hysteretic behavior, when an appropriate procedure is 
followed. It appears below a certain temperature Tc. 
We will propose that we have observed a transition from 
a vortex fluid stated to a new state occurring below 
Tc, which is characterized by the appearance of the hys- 
teretic behavior and of a critical velocity of the order ~ 1 
cm/s, beyond which the hysteretic component of NLRS 
is suppressed to zero. Actually hysteretic behavior itself 
has been reported by Kojima and his group at very low 
temperatures 12 71 and by Chan's group [28], as well as by 
Reppy's group[29]; however, none of them discussed the 
hysteresis in connection to the transition from the vortex 



3 




V^^ [^m/sec] 

FIG. 3: "hysteretic" AATLi^S'^^,- NCRIF as well as /^Q^l, 
vs logVac- From the linear extension of the logVac de- 
pendence we obtain a critical velocity, G^^IO mm/s to sup- 
press the NCRIF^pss to zero, which compares well with 
K:=h/(m4.^o - 7r)='^6-12 mm/s, for ^o=50-25 nm, see text. 



fiuid(VF) state. The VF state has been recently experi- 
mentally clarified by our report of the unique Vac and T 
dependent TO responses [9], effects which are absent in 
known superfluid transition responses. 

We describe here the first systematic study of the hys- 
teretic component of the solid ^He TO responses as a 
function of Vac as well as T below Tc and discuss a pos- 
sible order parameter of the supersolid(SS) state below 
Tc. In Fig. [2] the TO period change as a function of Vac 
is shown for various T. All the data points for each T 
down to 80 mK follow a single line as was the case for 
our data in a previous publication above Tp [9] . The data 
points below 75 mK actually follow two lines. These 
two lines were produced as follows: The first measure- 
ments were performed at ]/ac=maximum set at T above 
500 mK, as "equilibrium" condition measurements. Then 
they were prepared by cooling down the sample to the 
lowest T~ 48 mK, and then warmed up to the desired 
T < Tc and measurements were performed by changing 
Vac step wise downwards over a long enough time, typ- 
ically 12 hours or longer for the whole sweep one way, 
to allow for any relaxation at this constant T. Then the 
other series of measurements were prepared by reversing 
the excitation Vac change upwards to the measuring exci- 
tation step- wise after equilibrium is reached at each step 
up to the maximum Vac- After completing measurements 
at a T then another set of measurements at a different T 
was performed in a similar manner. We observe in Fig. [2] 
that the hysteretic component (the difference between the 
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FIG. 4: The T dependence of the NCRIF =/?ss(at 200 
/xm/s) (closed symbols) appeared as a result of the "hysteretic 
process" in % of the total He mass. It represents the "persis- 
tent" NCRIF produced by the process. It would be evidence 
of 3D macroscopic phase appearance below Tc. It compares 
with "paramagnetic" behavior of NLRS(Vac=0, open sym- 
bols), which follows the Langevin function as in Fig. 1, and 
suggests a "susceptibility" feature of NLRS. 



two passages of Vac) has a unique Vac dependence and in- 
creases towards lower T. 

Fig. [3] shows much more clearly the hysteretic compo- 
nents of both dissipation and NLRS as a function of Vac 
presented on a logarithmic horizontal scale. An inter- 
esting observation is that the hysteretic component ap- 
pears only above a characteristic AC velocity^ 40/im/s 
and it reaches a maximum at ~ 200/im/s. Furthermore, 
at higher AC velocity,^ 500/im/s it starts to decrease. 
It looks as if it follows a logVac linear relation passing 
through ANLRShys= in the range of Vac, 8i2 mm/s, 
signaling the total depression of the hysteretic NLRS 
component. The upper column of Fig. [3] shows the en- 
ergy dissipation change across the hysteresis loop AQJ^^^ 
as Vac A most impressive thing is that while AQ^J^ 
is negative ANLRShys is increasing over a considerable 
range of Vac- We suggest this is evidence for a shielding 
current preventing penetration of vortex lines as well as 
thermally excited vortices into the volume as in the Lan- 
dau state of superfiuids or Meissner state of superconduc- 
tors. At higher AC velocity above 500/im/s introduction 
of vortex lines would cause a logVac linear decrease of 
supersolid density pss- 

Fig. [4] indicates the T dependence of the hys- 
teretic component of AA/'Li^S'^^^ (200/im/s) together 
with "equilibrium" NLRS{Vac^O). The former appears 
below about 75 mK and increases almost linearly toward 
lower T to ~ 60 mK, then further increases more steeply 
as T lowers. We have tried to fit it with the expected be- 
havior for a 3D supersolid density, pss= (1 — T /Tc)^ =V 
with 7=2/3 and it yields an extrapolation to K(t=l). 
We can safely estimate " the supersolid density extrap- 
olated to K" within some error bar. It would be 
0.065±0.025% of the total solid He sample mass with 
Tc =56.7 mK. The pss{T = 0) is less than but on the 
same order as the expectation value of NLRS(T = 0), 
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of their early work on hysteresis in [3l| . 



FIG. 5: Critical behavior of ^ obtained from the data in Fig. 4. 
We took Tc — 56.7 mK, supposing pss— (1 — T/Tc)^—f^ with 
7=2/3. And we obtain ^o'^ 25 to 50 nm by simple extrapola- 
tions, horizontal and straight extension of the linear relation. 



which increases towards K following a Langevin func- 
tion, as discussed. From the absolute value of the pss^ 
we can evaluate super solid coherence length ^, following 
a common practice based on consideration of Josephson's 
length [3a]. We could fit rather weh with the ^-^/^ de- 
pendence expected for a 3D superfluid, obtaining an ex- 
trapolated ^ value at T=0, ^o=25 to 50nm as shown in 
Fig. m assuming 7714 to be the atomic ^He mass. We 
find yc=h/(m4.^o-7r)=~6-12 mm/s, for ^o=50-25 nm, a 
surprising coincidence. It is very interesting considering 
the microscopic origin of the 3D macroscopic supersolid 
phenomenon as suggested by the present work. Actually 
there have been various observations of the size of NLRS^ 
but nobody else discussed the transition from the VF to 
the SS state except the present authors. A similar hys- 
teresis phenomenon was already experimentally observed 
for a ^He sample with 30 ppm ^He long ago in an acoustic 
experiment by Iwasa and Suzuki in 1980 [31], but without 
noticing the relationship to a vortex state. 

We have found the start of the hysteretic behav- 
ior below Tc, and evaluated the hysteretic component 
ANLRShys cts a function of Vac and discussed it as su- 
persolid density Pss(T) and discussed the 3D coherence 
length ^ and a consistent critical AC velocity Vc^l cm/s. 
So far we have neglected the anisotropy of the hep crys- 
tal of "^He and treated it as an isotropic supersolid. The 
linear T dependence of pss {T) between 75 and ~ 60 mK 
may have some connection to this problem. 
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